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We propose a scheme that allows to coherently extract cold atoms from a reservoir in a determin- 
istic way. The transfer is achieved by means of radiation pulses coupling two atomic states which 
are object to different trapping conditions. A particular realization is proposed, where one state has 
zero magnetic moment and is confined by a dipole trap, whereas the other state with non-vanishing 
magnetic moment is confined by a steep microtrap potential. We show that in this setup a prede- 
termined number of atoms can be transferred from a reservoir, a Bose-Einstein condensate, into the 
collective quantum state of the steep trap with high efficiency in the parameter regime of present 
experiments. 

PACS numbers: 03.75.-b,03.75.Nt,32.80.Bx 



I. INTRODUCTION 

Recent experimental progress with ultracold atomic 
gases has opened exciting directions in the study of many 
body systems, aiming at the full coherent control of struc- 
tures of increasing complexity 0. Applications, like re- 
alizations of quantum information protocols, are very 
promising and presently investigated li iSs. iJ] • More- 
over, these systems allovif to study the frontiers between 
the quantum and the classical world ^ . 

In this perspective increasing interest has lately been 
attracted by the possibility of disposing single or few 
cold atoms by deterministic extraction from a reservoir, 
thereby allowing to control and manipulate them. This 
realization of tweezers applied to quantum objects has 
been called quantum tweezers. A recent proposal ex- 
ploited tunnelling from a condensate to a moving quan- 
tum dot in order to realize this scenario 

In this work we discuss an implementation of quantum 
tweezers, where trapping potentials, whose steepness de- 
pends on the atomic magnetic moment, are coupled by 
means of radiation. This procedure has been utilized for 
manipulating cold atomic clouds 0, ll| , and it has been 
proposed for engineering collective states of atoms in op- 
tical lattices 9]. Here, it is applied in order to coherently 
load atoms from a condensate into a steep trap by suit- 
ably coupling the atom electronic states. A sketch of 
the setup is shown in Fig. ^ and can be summarized as 
follows. Atoms are initially in a hyperfine state of zero 
magnetic moment, which we denote by \b), where the 
atomic center-of-mass is confined by a dipole trap, and 
form a Bose-Einstein condensate (BEG). The state \b) is 
coupled with radiation to the hyperfine state |a), which 
has non-vanishing magnetic moment and is subject to 
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FIG. 1: Extracting atoms on demand from a Bose-Einstein 
condensate (BEC): Atoms are transferred from the ground 
state of a BEC, confined by the potential Vb, into the one- or 
two-atom ground state of a steep potential Va (the tweezers). 
The horizontal lines indicate the corresponding energy levels 
(in arbitrary units) . The transfer is implemented by radiative 
coupling the electronic states |6) and |a), undergoing different 
trapping conditions. The number of atoms transferred into 
the tweezers is controlled by spectrally resolving the energy 
splitting between the one- and two-atom ground states of the 
tweezers. 



the steep potential of a magnetic trap. The atoms are 
selectively transferred by means of radiative coupling be- 
tween the condensate and a collective quantum state of 
the tweezers. The two traps can be found in the same 
spatial region, and then coherently separated spatially 
after the transfer has occurred |a.llO|. 

We investigate the dynamics and efficiency of quantum 
tweezers using radiation and show that their efficiency for 
transferring a certain number of atoms into the quantum 
state of a trap can be larger than 99% in experimen- 
tally accessible parameter regimes mill El- We study 
three implementations of atom transfer into the quan- 
tum tweezers: (i) by means of the adiabatic passage by 
adiabatically ramping the radiation frequency across the 
resonance |l4l | (ii) by combining adiabatic and diabatic 
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passages by means of laser pulses ^^'^ (iii) t>y res- 
onant coupling using pulses of well defined area. The 
various techniques are compared and their efficiency is 
discussed in the parameter regime of experiments with 
microtraps [l6j . 

This article is organized as follows. In section ITU the 
theoretical model is introduced and the relevant param- 
eters are identified. In section IIIII the efficiency of the 
transfer for different realizations of the radiation cou- 
pling between the two hyperfine states is discussed. In 
section IIVI the conclusions are drawn and outlooks are 
discussed. In the appendix a detailled analysis of the 
parameter regime, where the quantum tweezers can be 
implemented, is presented. 



II. THEORETICAL DESCRIPTION 

We consider N ultracold atoms of mass M whose rele- 
vant internal degrees of freedom are the stable hyperfine 
states I a) and \b). The atomic center of mass experiences 
at position x = (x,y,z) a harmonic potential Vj{x) of 
the form 
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(1) 



where Vjx^ i^jy, are the frequencies along each carte- 
sian direction and j — a,b labels the atomic state. The 
potential 14 (x) is steeper than Vh(x), namely the fre- 
quencies Vax,Vay, Vaz > miu ( z/fea; , j/fej^ , V},^). At time t = 
the atoms are all in the state \b) and constitute a Bose- 
Einstein condensate. The aim is to coherently transfer 
a predetermined number of atoms from the condensate, 
which acts as a reservoir, into a collective quantum state 
of the atoms in the steep potential Va (x) . 

Transfer of atoms into the quantum tweezers are im- 
plemented by suitably driving the states \a) and \h) with 
radiation pulses, which can be in the microwave or op- 
tical regime, thereby implementing a Raman transition. 
Coherent and selective transfer is achieved provided the 
radiation pulses are sufficiently short, so that incoherent 
dynamics during the interaction is negligible, and suffi- 
ciently long to spectrally resolve a selected atomic transi- 
tion. This latter statement will be quantified in Sec. Ill Bl 
In the frame rotating at the atomic dipole resonance fre- 
quency, the relevant dynamics describing the interaction 
of the atom with radiation is given by the Hamiltonian 
7i, which we write as 



H = Tia +'Hb + Ti-sc + 'Hint{t) ■ 



(2) 



Here, the terms Ha and Tib describe the dynamics of the 
atom's center-of-mass in the states \a) and respec- 
tively, and have the form 



^2 

v2 + y,(x) 



9jj / c^x?/;j(x)'f7/;j(x)'''-0j(x)V'j(x) , (3) 



where ('0j(x)^) is the field operator annihilating 

(creating) an atom at position x in the internal state 
\j), the term gjj = Airh^ajj /M represents the interaction 
strength of two-body collisions of atoms in the state \j) 
and ajj denotes the s-wave scattering length (j = a, 6). 
Collisions between atoms in different hyperfine states are 
described by the term 

Use = gab / dxi/)a(x)V6(x)'^?/'a(x)V'6(x) (4) 



where gab = '^T^h^cLab/M is the interaction strength, with 
the s~wave scattering length aab for collisions between 
one atom in state \a) and one atom in state \b). The cou- 
pling with radiation is described by the time-dependent 
term 

Wi„t(i) = -ftA(f) / dx^,(x)t^,(x) 



+ -h^lLit) j dx (V'a(x)V6(x) -I- B.C.) , (5) 

where ^L{t) is the real- valued Rabi frequency and A(i) is 
the detuning of the radiation at time t from the resonance 
frequency of the transition \b) — s- \a). 

All atoms are initially in the internal state |6) and form 
a condensate. We denote this collective state, the ground 
state of Tib for N atoms, by |0). By means of a radiation 
pulse a number n of atoms (with n <^ N) is transferred 
into the lowest energy state of Ha- We denote by \n) the 
state of the system after the transfer, corresponding to 
N — n atoms in the ground state of Tib and n atoms in 
the ground state of Tia ■ The efficiency of the procedure is 
measured by the probability Pg^niT), which corresponds 
to the probability at time T, at the end of the pulse, of 
finding n atoms in the ground state of Tia- It is defined 



with 



(6) 



(7) 



Po^„(T) = \{n\U{T)\0)\' 



U{T) = T exp (j^ dTH(T)^ 



and T indicates the time ordering. In section IIIII we 
describe strategies of varying the coefficients A(t) and 
^L{t) in Tiint as a function of time in order to achieve 
unit efficiency. 



A. Basic assumptions and approximations 

In order to study the dynamics and the efficiency of the 
tweezers, we decompose the field operator in a convenient 
basis. Here, we use the Fock decomposition of the field 
operator ■(/'a(x), namely 



V'a(x) 



E 



(8) 
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where fia labels the excitations of the harmonic oscillator 
with Hamiltonianp^/2M+Va(x), ^JJ" (x) is its wave func- 
tion and afi^ is the field operator creating an atom in the 
corresponding state |r?a). We remark that the state {fta) 
describes single particle eigenstates, hence they are not 
the eigenstates of Tia when two or more (interacting) 
atoms are inside the tweezers. Nevertheless, in this case 
it still provides a plausible ansatz for estimating the order 
of magnitude of the physical parameters. 

Since the number of atoms which are extracted from 
the condensate is negligibly small compared with N, 
we replace the field operator for atoms in the hyperfine 
state \b) by a scalar. At t = the atoms are assumed 
to be in the condensate ground state with macroscopic 
wave function '/>6(x) and chemical potential fi. These 
quantities are assumed to be constant and negligibly af- 
fected by the transfer of atoms into the tweezers. The 
condensate is hence a reservoir of atoms at given den- 
sity nb(x) = A^|(/if,(x)|^. The condensate excitations are 
accounted for by the (low-energy) collective modes at fre- 
quencies Wg, and the part of Hamiltonian Tii, which may 
limit the efhciency of the quantum tweezers takes the 
form Hb = ^q^i^qbqbq where bq, 6 J are the annihila- 
tion and creation operators of a phonon of energy hujq. 
In addition, quantum tweezers overlap spatially with the 
condensate, and collisions between the atom inside the 
tweezers and the condensate, described by Eq. can 
be detrimental for the coherence of the process. Their in- 
fluence on the tweezers efficiency is discussed in Sec. Ill Bl 
Below, we treat the Hamiltonian term Q as a small per- 
turbation, whose effect is to introduce small shifts of the 
energy levels. 

In the situations we consider the condensate ground 
state is coupled (quasi-) resonantly with the ground state 
of Tia for n-atoms and one or two atoms at a time are 
transferred into the steep potential. The states which are 
relevant to the dynamics are denoted by |0), |1) and |2), 
corresponding to 0,1,2 atoms in the ground state of the 
quantum tweezers, respectively. Their energy is E{n) — 
{n\H\n) with n = 0, 1, 2. By setting £'(0) = 0, they take 
the form 



and 



E{n) = Ea{n) + Esc{n) - nfj,, 



(9) 



where Ea{n) is the energy of n atoms in the ground state 
of the Hamiltonian TCa and Esc{n) is the interaction en- 
ergy due to the term while the term — n/i accounts 
for the extraction of n atoms from the condensate. The 
explicit expressions for Ea{n) and Esc{n) are found from 
the overlap integral between the condensate wave func- 
tion and the wave function describing the ground state of 
Tia for n atoms. We evaluate them by approximating the 
latter with the wave function of n non-interacting atoms 
in the ground state of the oscillator, and find 



Eain) = n 



\ j=x,y,z 



hA 



(10) 



Esc{n) g^f,nN / dx |(?!)a(x)p|(/)f,(x) 



(11) 



Here, we have denoted by (/)a(x) the state (/>"°(x) with 
fia = (0, 0, 0) and AEn is an energy shift due to particle- 
particle interaction in state |a}, which can be decomposed 
into the terms AEn ~ A_Ecoii + AEdipo\c{t)- The term 
AE'coii accounts for the collisions of atoms in the ground 
state, which we estimate for n = 2 as 



A-Bcoll ~ gaa j dx \4>a{y^)\ 



(12) 



while the term A£'dipoie(0 due to dipole-dipole radia- 
tive interaction between the atoms inside the tweezers 
and it is proportional to the dissipative rate associated 
with the radiative transfer [T8l| . 

The coupling due to radiation between the states \n) 
and |n -f 1) is described by the matrix elements (n|7i|rt-|- 
1) ^ [n + l\n\n)\ with n = 0, 1. We denote by 



f7W = -(n|7^|n+l) 



(13) 



the corresponding Rabi frequency, which takes the ex- 
plicit form 



rfl^) = nL^N{n + l) / dx0fe*(x)0„(x) . (14) 



We assume that atoms are extracted from the center of 
the condensate, which here corresponds with the center 
of the tweezers. For a very steep tweezers trap, whose 
size is much smaller than the size of the condensate, the 
integral in Eq. (|14|l is approximated by 



17(1) ^ nL^nb{n+l) / dx0,(x) 



(15) 



where Ub — nb{Q) is the density at the centre of the con- 
densate. 

The resonant coupling between the states |0) and |2) 
can be characterized by the Rabi frequency 



(2) 



E2/2h ' 



(16) 



where the denominator corresponds to the value of the 
detuning for which the two states are resonantly coupled, 
and 

E2^ f^'^aj + ^E2 + E,,{2) - 2fi . (17) 

j=x,y,z 

The coupling, Eq. H16() . has been evaluated in the limit 
ri^^) <^ E2/h in second order perturbation theory p^ . 

For later convenience, we denote by A = Ei/h the 
detuning for which the states |0) and |1) are resonantly 
coupled, with 



(18) 
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B. Values of the parameters 

In this section we estimate the order of magnitude of 
the energy terms and check the consistency of our as- 
sumptions for experiments where the tweezers trap is a 
magnetic microtrap. In particular, we consider the pa- 
rameter regimes discussed in |l6j| . 

We first focus on the spectrum of one atom inside 
the tweezers. According to Eqs. (Pl- ljTU)) . the energy of 
the vibrational ground state |1) reads E{1) = Ea{l) + 
■Esc(l) - M with 



(19) 



For Va — liTX 30 kHz and vi, = 2tix 100 Hz and a conden- 
sate of TV w 10^ *'^Rb atoms with density rih sa 3 • 10^^ 
atoms/cm'^, one finds fi/h k, 1.8kHz while i?sc(l)/ft is 
orders of magnitude smaller. For the given value of Ua 
the size of the tweezers ground state wave function is 
flTweezor ~ 60 nm. Morcovcr, the energy distance be- 
tween the ground and first excited state is of the order of 
hva- Thus, the frequency of the steep trap sets a funda- 
mental limit for selectively addressing the ground state of 
the one-atom tweezers, thereby avoiding the excitation 
of other single-atom states. 

The ground state of two atoms inside the tweezers has 
energy E{2) = Ea{2) + E^^{2) - 2/i with 



Ea{2) = ihva - 2hA + AE2 



(20) 



where £;sc(2) 2£;sc(l). The term AE2 = AE^coii + 
Ai?dipoio describes the energy shift due to atom-atom 
interaction inside the tweezer. The term due to collisions 
is AiJcoU /h^ 27r X 2 kHz. When the coupling between the 
hyperfine states is achieved with coherent Raman transi- 
tions and radiation in the optical range, the term due to 
dipole-dipole interaction is of the order of the excitation 
rate of the state \a) and it is negligible for the consid- 
ered setup. Hence, the frequency of the steep trap is the 
largest frequency scale, which determines also the size 
of the energy separation between the ground and first 
excited state of the two-atom tweezers. 

From Eq. (|10|) it is visible that selective addressing of 
the states |1) and |2) requires to resolve frequencies of 
the order of AE2/h, namely the interaction energy of 
the atoms inside the tweezers. This sets a fundamen- 
tal limit on the parameters for implementing the trans- 
fer, and consequentely on the time for implementing the 
transition. 

As in this parameter regime AE2 S> -Esc, it is justified 
to treat the Hamiltonian term, Eq. (0J, as a small pertur- 
bation. Neglecting the coupling of the state \a) with the 
excitations of the condensate is more delicate. This cou- 
pling can be neglected when the pulses spectrally resolve 
the condensate excitations. This assumption may be rea- 
sonable when the trap frequency > 2tt x I kHz. In this 
work we will neglect this coupling, assuming a parameter 
regime where these excitations can be spectrally resolved. 



For smaller values of Vb, however, its effect must be taken 
into account. Its systematic study in such regime will be 
object of future investigations. 

We discuss now the approximation on the ground state 
wave function of |n). It is in fact questionable whether 
one can approximate the wave function of the state 
with n > 2, with the product of single particle wave 
functions. We note, however, that the size of the ground 
state wave function in the steep traps we consider here 
is still larger than the s-wave scattering length. Numer- 
ical studies have shown that two-body s-wave scatter- 
ing is still a reasonable approximation in relatively steep 
traps [23, 1^ I23I . Hence, although the ansatz we use is 
not reliable for exact results, it is still plausible for gain- 
ing insight into the efficiency of the tweezers. Neverthe- 
less, it is understood that the application of the schemes 
discussed below to a certain experiment requires the ac- 
curate knowledge of the relevant parameters, which have 
to be evaluated for the specific atomic species and trap- 
ping conditions. 



III. TWEEZING ATOMS WITH LASERS 

In this section we discuss three methods of shaping 
radiation pulses in order to transfer a definite number of 
atoms from the condensate into the microtrap. The first 
method uses adiabatic ramping of the frequency of the 
radiation coupling |rt — 1) — > \n) (adiabatic passage) 
14]. The second method is based on the combination 
of two pulses, which induce a sequence of adiabatic and 
diabatic passages (SCRAP) The third method uses 
resonant pulses with well defined pulse area [l^. 



A. Adiabatic Passage 

We consider the transfer of n atoms from the conden- 
sate into the quantum tweezers by adiabatically ramping 
the frequency of radiation coupling the two atomic states. 
The frequency of the pulse, namely the detuning A(t) in 
Eq. 0, is slowly varied as a function of time from the 
initial value A.^ to the final detuning A f , allowing the 
system to follow adiabatically the evolution. 

Figure|21displays the energy eigenvalues of W as a func- 
tion of the detuning A. Adiabatically sweeping the de- 
tuning through these resonances corresponds to remain 
in an instantaneous eigenstate of the Hamiltonian Ti.. 
Transfer of one atom is implemented by ramping the 
detuning between the values A^ and A/ in Fig. |21 cor- 
responding to change the energy eigenstate along the 
dashed curve. Adiabatic evolution is preserved when the 
frequency ramping rate A fulfills the relation 



|A(to)| « 



n{AE2/hy 
2|ln(l-Po) 



(21) 



where AE2 is the two-atom interaction energy and Pq is 
the threshold transfer probability, such that for Pq^i > 
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FIG. 2: Extraction by adiabatic ramping the laser detuning A 
through atomic resonance. The color lines show the energies 
E{n), when n atoms are inside the tweezers ground state, as a 
fimction of the detuning. The energies E(n) cross the energy 
E{Q), horizontal line, at different values of A, corresponding 
to the interaction energy of n atoms inside the tweezers. The 
dashed lines are the dressed state energies, obtained by radia- 
tive coupling between the BEG and the tweezers. When the 
atoms are all initially in the BEG and the detuning is ramped 
adiabatically, the system dynamics follow correspondingly the 
lower dashed line. Transfer of one atom into the ground state 
of the tweezers is achieved by ramping through the avoided 
crossing between E{Q) and E{1). 



Pq the transfer into the tweezers has been successfully 
implemented. The derivation of inequaUty H21|) is shown 
in Appendix A. 

The upper quadrants of Fig. |31 display the transfer ef- 
ficiency Pq— »i as a function of the rate of ramping | A| for 
two values of the microtrap frequencies = 27r x 30 kHz 
and Va = 2tt X 100 kHz. Here the red curve has been ob- 
tained from numerical simulations, and the shaded region 
denotes the range where the transfer probability is larger 
than 0.99. The black curve corresponds to the prediction 
of the Landau-Zener formula in the adiabatic reg ime, 
namely for large transfer efficiencies [23 . [25l |26| and 
discussed in appendix A. This takes the form 



(LZ) _ 



1 



(22) 



with aad = 7r5^/2?i^|A(io)| the adiabaticity parameter 
and A (to), 5 the rate of ramping and the energy split- 
ting, respectively, at the avoided crossing between E{Q) 
and E{\). Efficiencies of the order of 99% are achieved 
for transfer times of the order of the millisecond. The re- 
sults show that faster transfer is achieved at larger trap 
frequencies. In fact, the interaction energy /S.E2 increases 
with Va, thereby allowing for larger values of the ramping 
rate according to Eq. H21|l . 

The transfer of two atoms into the tweezers can be 
achieved starting from the state |0) by adiabatically 
ramping the detuning A either (i) subsequently across 
the resonances £'(0) = £'(1) and then E{1) = E{2) 



FIG. 3: Transfer efficiency Po->i (top) and Po^2 (bottom) as 
a function of the ramping rate |Aj for (a) Va = 2% x 30 kHz 
and = 4 kHz (b) Va = 2-KX 100 kHz and f^L = 30 kHz. The 
red curves show the numerical result, the black curves display 
the theoretical predictions of the Landau-Zener formula. The 
shaded areas denote the regions, where the theoretical success 
probability, according to Eq. 12211 . is above 99%. Efficient 
transfer can be achieved in few milliseconds. 



(lower dashed line in Fig.[2Jl, or (ii) across the resonance 
E{Q) — E{2). The first case corresponds to transfer- 
ring one atom at a time. In the second case two atoms 
are transferred simultaneously. This latter implementa- 
tion has a drawback due to the small value of the radia- 
tive coupling between the states |0) and |2). In fact, the 
corresponding Rabi frequency, estimated in Eq. (|16|l . is 

in general rather smaller than the Rabi frequency 
Eq. (|14|) . describing the transfer of one atom into the 
tweezers. This smaller values leads to smaller energy 
splitting at the avoided crossing between the levels |0) 
and 1 2), therefore imposing much stricter limits on the 
ramping speed. The lower quadrants of Fig. |21 depict 
the results of the simulation for the sequential adiabatic 
transfer of two atoms into the microtrap. The efficien- 
cies and transfer times are comparable with the effi- 
ciencies reached for transferring a single atom, showing 
that two atoms can be transferred into the ground state 
of a tweezers trap with ^'^ = 27r x 30 kHz in few mil- 
liseconds. The black curve corresponds to the Landau- 
Zener prediction, which in the adiabatic region is to 
good approximation the product of the transfer probabil- 
ities at the single avoided crossings, and takes the form 

(LZ) _ p(LZ) p{LZ) jj^^^ p^^^ j^^^^ ^^^^ 



^0^2 ~ o->i ^1^2 



where the energy splitting 6 is h 



B. Adiabatic Passage using Pulses 

We discuss now transfer of atoms into the tweezers by 
means of the so-called Stark-chirped rapid adiabatic pas- 
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FIG. 4: Sketch of the dynamics when implementing the 
SCRAP technique for transferring one atom into the tweez- 
ers, (a) Pump pulse (f2_L(t)) and a.c.-Stark shift (Ai(t)) as 
a function of time, (b) Corresponding spectrum. The color 
lines in (b) are the energies E{n) as evaluated at the cor- 
responding value of Ai(t). The dashed curves indicate the 
dressed state energies. The vertical solid (dashed) lines indi- 
cate the instants at which the transition |0) |1) (|0) |2)) 
is driven resonantly. 



sage (SCRAP) technique [2^. Here, population transfer 
between two quantum states is achieved by suitably vary- 
ing the Rabi frequency ^Lit) and the detuning A(t) in 
Eq. jSJ, thereby combining adiabatic and diabatic tran- 
sitions. In this scheme, two pulses of finite duration are 
utilized, which are time-delayed one with respect to the 
other. One pulse is far-off resonance and induces an a.c.- 
Stark shift A(<) on the transition \b) \a), while the 
second pulse, which we denote by 'pump pulse', couples 
(quasi-)resonantly the states \a) and |6) and has Rabi 
frequency 0/, (t) . The detuning is shaped in order to ful- 
fill the resonance condition A{t) = Ei/h at two instants 
of time ti and t2 ■ The shaping and time-delay between 
the two pulses is such that at ti the system undergoes an 
adiabatic passage of the same type as the one discussed in 
Sec. IIII Al while at t2 the transition is diabatic, thereby 
ensuring that at the end of the pulses the transfer be- 
tween the states |0) and |1) is achieved |2^. A possible 
realization of the pulses as a function of time and the cor- 
responding dynamics are displayed in Fig. 0] Here, the 
energy levels E{Q) and E{1) cross at the instant ti — 
and t2 — T. At these instants, the system undergoes 
an adiabatic and diabatic transition, respectively. As a 
result, if initially all atoms are in the condensate, then 
finally one atom is transferred into the ground state of 
the tweezers. 



Since at least at one avoided crossing the dynamics 
must be adiabatic, the bounds on the values of the Rabi 
frequency and on the time variation of the a.c.-Stark 
shift are basically the same as for the adiabatic passage 
in Sec. IIII Al Moreover, the finite size of the pulses and 
the requirement that a diabatic transition takes place 
at one of the crossing, introduces a further parameter, 
namely the width of the pump pulse Tji, on which the 
transfer efficiency depends. The parameter regime for 
the applicability of this technique is discussed in detail 
in Appendix B. 

Figure [Sja) displays the transfer efhciency as a func- 
tion of Tfj and Q^, the maximum Rabi frequency of the 
pump pulse, as obtained from a numerical simulation for 
a tweezer trap with i/q = 27r x 30 kHz. Transfer effi- 
ciencies Pq^i > 99% are achieved in the white region, 
corresponding to transfer times of the order of some mil- 
liseconds. High efficiencies are reached for a relatively 
wide range of parameters, showing that the method is 
robust to fluctuations of the parameters of the pulses. 
These datas have been evaluated under the condition for 
which the pump pulse achieves its maximum value at 
the same instant of time at which the detuning pulse 
fulfills the resonance condition, as depicted in Fig. 2| In 
this case the ideal conditions for adiabaticity are reached. 
The effect of a time-lag 6t between these two events is 
shown in Fig.l^Jb), which reports the dependence of the 
transfer efficiency on St. The transfer efficiency exhibits 
a wide plateau around St = 0, which is of the order of a 
millisecond, showing that the method is also very robust 
against this kind of fluctuations. Note that the curve ex- 
hibits a second maximum at St « —4 ms, corresponding 
to the case in which the pump pulse is centered at the 
second instant t2 at which A fulfills the resonance con- 
dition A = Ei/h, namely to the situation in which the 
sequence of diabatic and adiabatic passages is reversed. 
Ideally, the two maxima should be symmetric. Asymme- 
try here arises from the fact that the dynamics at the 
crossing corresponding to the resonance for simultaneous 
transfer of two atoms into the tweezers (dashed verti- 
cal line in Fig. Efb)) are not perfectly diabatic. This 
introduces a fundamental difference in the efficiency of 
transfer, which depends on whether this transition point 
is crossed, namely on whether the diabatic transitions 
occur before the adiabatic one. 

Similarly to the adiabatic passage in Sec. IIII Al trans- 
fer of two atoms by means of SCRAP is better imple- 
mented sequentially, i.e., by transferring one atom at a 
time. Figure shows a pulse sequence and the corre- 
sponding energy spectrum as a function of time for trans- 
ferring sequentially two atoms inside the tweezers. The 
ideal dynamics follow the lower dashed line in Fig.|Hfb). 
The transfer efficiency as a function of Tq and Q is dis- 
played in Fig. [7| Here, efficient transfer of two atoms, 
larger than 99%, is achieved in times of the order of sev- 
eral milliseconds. The relatively broad range of values 
for which Po^2 exceeds 99% shows that the method is 
robust against parameter fluctuations. 



0.0 0.5 1.0 1.5 2.0 

Th/ms 
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-4.0 -2.0 0.0 2.0 

St I nis 



FIG. 5: (a) Contour plot for the transfer probability Po-,i as 
a function of the coupling pulse maximum Q. and width Tn. 
The tweezers trap frequency is = Ztt x 30 kHz and the a.c- 
Stark shift pulse has width Ta = ITq,. Pump and Detuning 
pulses have Gaussian envelope. The red cross indicates the 
parameters Vl = 15kHz and Tn = 1ms, which are used in 
evaluating the curve in (b). (b) Probability Po^i as a function 
of the time variation &t between the instant of time at which 
the resonance condition is achieved by the a.c.-Stark pulse 
and the instant of time at which the pump pulse reaches its 
maximum value. The maximum at &t « —4 ms corresponds 
to the dynamics for which the temporal sequence of adiabatic 
and diabatic passages is exchanged. The asymmetry between 
the two peaks is discussed in the text. 



C. Transfer by population inversion 

We finally discuss transfer of atoms into the tweez- 
ers by radiation pulses of chosen areas, coupling reso- 
nantly the states |0) and Perfect transfer is achieved 



1 1 ms 

FIG. 6: Sketch of the dynamics when implementing the 
SCRAP technique for transferring sequentially two atoms into 
the tweezers, (a) Pump pulse (f2i,(t)) and a.c.-Stark shift 
(Ai(i)) as a function of time, (b) Corresponding spectrum. 
The color lines in (b) are the energies i5(n, t) as evaluated 
at the corresponding value of Ai(t). The dashed curves in- 
dicate the dressed state energies. The vertical lines indicate 
the instants at which the transition |0) jl) and |1) |2) 
are driven resonantly. Transfer of two atoms is achieved when 
the dynamics start at the blue solid curve. 



when /S. = Ei/h and when the the pulse area fulfills the 

relation / dr Tl'^i^ (t) = tt . The limitations are on the 
choice of the Rabi frequency, whose value is bound by the 
interaction energy in order to have negligible coupling to 
off-resonant states (see Appendix A), and on the pulse 
duration, which must be sufficiently long to guarantee 
the spectral resolution of the tweezers energy levels. In 
Fig. |S| the transfer efficiency Pq^i is plotted for Gaus- 
sian pulses as a function of the pulse maximum intensity 
and duration. Efficiencies close to unity are realized with 
pulses of millisecond duration. 

Efficient transfer of two atoms is achieved by sequen- 
tially sending two pulses, each transferring one atom 
into the tweezers. Here, the first pulse is resonant with 
the transition |0) — > |1), thereby implementing a tt- 
rotation. A 7r-rotation on the transition |1) — > |2) is 
achieved by coupling it with a second pulse with detun- 
ing A — {E2 — Ei)/h and suitable pulse area. The spec- 
tral resolution, at the level of AE2/h, avoids that during 
the second pulse the atom in the tweezers is transferred 
back to the condensate, as the transition |0) — > |1) is 
far-off resonance. Transfer efficiency exceeding 99% are 
achieved on time scales of the order of several millisec- 
onds. 
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FIG. 7: Contour plot for the transfer probability Po->2 as 
a function of the coupling pulse maximum Q, and width Tq. 
The tweezers trap frequency is Va = '2n x 30 kHz and the 
a.c.-Stark shift pulse has width Ta ~ 2Tq. The detuning 
pulse has Gaussian envelope, the pump pulse is given by 



Ik 

2 



tanh(- 



tanh(- 



where tramp is the 



time in which the pulse is ramped up to the maximum value. 



IV. DISCUSSION AND CONCLUSIONS 

We have investigated the feasibihty of a scheme that 
employs radiation, thereby achieving the deterministic 
extraction of atoms from a Bose-Einstein condensate and 
their transfer into the quantum state of a steep trap. We 
have shown that a high transfer efficiency can be achieved 
in times of the order of some milliseconds in the param- 
eter regime of present experiments vifith microtraps |l6j | . 
The radiation used for implementing the extraction can 
be either in the microwave regime, thus coupling directly, 
say, a magnetic dipole transition like in |3,|8|, or in the op- 
tical domain, thereby implementing Raman transitions. 
Their application depends on the specific details of the 
system. 

The results here presented refer to ideal conditions, 
where at t = all atoms are in the condensate ground 
state and the coupling with the other excited states of 
the condensate can be neglected. The presence of non- 
condensed atoms and the coupling of the state \a) to the 
condensate excitations have not been considered. The 
presence of non-condensed atoms introduces a source of 
noise into the process, giving an indetermination in the 
transfer efficiency of the order of the non-condensed frac- 
tion. The coupling of the state \a) to the condensate 
excitations constitutes a further source of noise which 
limits the efficiency of the tweezer. Its effect can be 
eliminated by resolving spectrally the condensate exci- 
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FIG. 8: Transfer of one atom by resonant tt rotation, (a) Con- 
tour plot of Po-»i by means of a Gaussian pulses with A — 
E\/h as a function of the maximum pulse value 57 and of 
the pulse width Tn for i^a = Stt x 30 kHz. The black regions 
correspond to the parameter regime, where the efficiency is 
below 80%. (b) Po^i as a function of the pulse strength Q. 
for constant pulse width Tn = 1.5 ms. The corresponding 
parameters are indicated by the vertical line in (a). 



tations. This is possible when the condensate is confined 
in a sufficiently steep trap. In this case, the condensate 
excitations must be taken into account when setting the 
parameters, as they determine the spectral resolution to 
be achieved and eventually the transfer duration. 

Our analysis has been restricted to the transfer of one 
and two atoms into the corresponding tweezers ground 
state. Similar considerations apply for the transfer into 
one- and two-atom excited states. In this case, the pa- 
rameters must be carefully set, so to spectrally resolve 
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the desired level. The process, in this case, may take ad- 
vantage of symmetries, ruling out the coupling to states 
which may be close but orthogonal to the initial state. 
The procedure discussed in this paper can be as well ex- 
tended to the transfer of three or more atoms into the 
tweezers. In this case, we expect that the scheme is effi- 
cient when transferring sequentially one atom at a time. 
When confining three or more atoms into the steep trap, 
however, one must consider many-body effects which can 
give rise to instabilities. The dynamics and time-scales 
on which they manifest themselves depend on the trap 
parameters and atomic species and has to be analysed 
case by case. 

To conclude, the possibility of extracting atoms on de- 
mand from a condensate realizes the deterministic cou- 
pling of an atom trap to a reservoir, thereby finding sev- 
eral analogies in quantum optical systems, and opens in- 
teresting perspectives in the study and the control of co- 
herent matter waves. 
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APPENDIX A: LIMITS FOR THE APPLICATION 
OF THE ADIABATIC PASSAGE 

In this appendix we discuss the parameter regimes for 
which the dynamics described in Sec. IIII Al apply. We 
first summarize the dynamics of the adiabatic passage, 
by discussing the transfer of one atom into the tweezers. 
The considerations made here can be directly extended 
to the transfer of n atoms. We assume that the detuning 
is swept between the initial value and the final value 
A/ in the interval of time T, such that A(0) = A^ and 
A(r) = A/. Efficient transfer from |0) to |I) is achieved 
when at these values of the detuning the states |0), |1) 
are eigenstates of Ti to good approximation. The value of 
the Rabi frequency is bound from the constraint that the 
relevant dynamics must involve only the states |0} and 
1). When these conditions are fuffiUed, the dynamics 
can be restricted to the subspace of states {|0) , |I)} and 
the eigenvalues of the reduced Hamiltonian at a given 
instant t, < t < T, are 



e±(I,i) = i/iAi(t)±i^/Ai(t)2 



(1) 



(AI) 



with 

Ai{t) = A{t) - Ei/h. (A2) 

The corresponding eigenstates take the form 

!$+(<)) = sine(t) |0) +cose(t) 1 1) (A3a) 
!$_(<)) = cos0(i) |0) - sine(t) 1 1) (A3b) 



where the mixing angle O is defined by the relation 



tan2e(t) 



(1) 



< e{t) < TT. 



(A4) 



We denote by to the instant of time at which the tran- 
sition is resonantly driven, namely when the detuning 
fuffills the relation A{to) = Ei/h. The minimum value 6 
of the energy splitting is reached at t = to and takes the 
form 



6=\e+il,to)^e.il,to)\^h n\ 



(A5) 



We now identify the parameter regime for which one may 
restrict the dynamics to the levels |0) and |I). This 
approximation is justified when (i) coupling to single- 
particle excitations of the one-atom tweezer is negligible, 

i.e. fl'^^ <C Va and (ii) when the state |2) contributes 

negligibly to the dynamics, namely when the condition 



(1) 



<C AE2/h 



is fulfilled. This condition on il^^-* ensures also a complete 
adiabatic transfer by ramping the detuning across the 
resonance E{0) = E{1) while keeping H.^ constant. In 
fact, in this limit the final detuning A/ can be chosen 
sufficiently far away from the resonance E2/2h, and at 
the same time be sufficiently large so that at A/ the 
unperturbed states |0) and |1) are to good approximation 
eigenstates of H. 

We now evaluate the probability Po^i of transferring 
one atom into the ground state of the tweezer when the 
dynamics can be restricted to the levels |0) and In 
the adiabatic regime the probability Po^i can be evalu- 
ated by using the Landau-Zener formula, which at the 
asymptotics, i.e., for large transfer efficiencies, takes the 
form 



(A6) 



1 



Here aad is the adiabaticity parameter, defined as 



t5^ 



dad 



2n^\A{to)\ ' 



(A7) 



(A8) 



The parameter ctad is proportional to the energy split- 
ting 6, defined in Eq. (|A5|I . and is inversely proportional 
to the time derivative of the detuning A at the level- 
crossing. Hence, the probability Pq^i approaches unity 
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for large values of Oad, i-e. for large values of the energy 
splitting and/or small values of the rate A. 

We fix now the threshold value Pq, such that for 
Pq^i > Pq the transfer is considered successful, and 
we evaluate the rate A, and thus the minimum time, 
required for transferring successfully one atom into the 
tweezers. From Eq. (I22|l we find that Pq^i > Pa when 



|A(io) 



< 



n{d/hy 



2|ln(l-Po) 



(A9) 



Using Eqs. IjASp and 



61 we rewrite this condition as 

2 



\A{to)\ < 



7r{AE2/hy 
2|ln(l-Fo) 



(AlO) 



Hence, the time ri needed for transferring one atom into 
the tweezers must fulfill the relation 



chosen ti < t2, so that the adiabatic passage occurs at 
the first crossing, and set ti = and t2 = t. 

The system dynamics can be reduced to the two levels 
0) and |1) provided that the relation in Eq. (|A6|I is ful- 
filled, which set an upper bound to the maximum value 
fl that the pump pulse reaches. The dynamics are adia- 
batic around the point ti provided that the pulses time 
variation is sufficiently smooth such that the adiabaticity 
parameter in Eq. (|A8|) is very large. For this reason, ide- 
ally the pump pulse should reach its maximum at <i, such 
that the level splitting is maximum at the avoided cross- 
ing. The adiabaticity of the transfer at ti can be eval- 
uated by checking for which parameters the adiabaticity 
parameter a^d of Eq. (|A8(I is much larger than unity. 
From Eqs. (|A6|) and IjASp one finds a condition for the 
temporal variation of the detuning. For Gauss pulses as 
in Eq. (|Blb|) this condition takes the form 



n = 



A(io) 



> 



In(l-Po) 



AE2/h 



(All) 



where we have taken |A/ — Ai| AE2/h. With the 
parameters of section Hi Bl for AE2/h = 27r x 2 kHz and 
a success probability Pq = 99%, the transfer time must 
fulfill the relation ti ^ 0.2 ms. 



APPENDIX B: LIMITS FOR THE APPLICATION 
OF SCRAP 

In this appendix we discuss the parameter regimes for 
which the dynamics described in Sec. IIII Bl applies. For 
simplicity we consider the case in which both the pump 
as well as the far-off resonance pulse have Gaussian en- 
velope, and are given by 



ilL{t) — fiexp 



2 1 



and 



^offset 



t 



A cxp 



(Bla) 



(t-rf 



(Bib) 



where Tq and Ta are the pulse widths, the terms and 
A denote the pulse maxima, and the offset value AoSsct 
is such that at the instants ti and t2 the resonance condi- 
tions A{ti) — A{t2) — Ei/h are fulfilled. Here, we have 



nAEl 
4 ft2 



At 

^ A 



rp2 



(B2) 



where here to — ti — 0. 

Ideally, the adiabatic transition between two eigen- 
states of Ha + Ti-b, say |0) to the state |1), is imple- 
mented over an infinite time. As the pump pulse has 
a finite duration Tfj, then To must exceed a minimum 
value tjump in order to achieve sufficiently large transi- 
tion probabilities at the instant ti. Using the relation 
^jurnp 2rii/|A(ii)| [23|, where fli is the largest value 
of the coupling between |0) and |1), this condition corre- 
sponds to 



Tc 



>f^(i)Z|exp 
tA 



rp2 



(B3) 



which fixes a relation between the widths and maxima of 
the two pulses and their relative delay. Simple algebra 
shows that this condition is consistent with the condition 
in Eq. provided that AE2Tn/h > 1. 

Finally, the transition is diabatic aX t2 ^ t when the 
pump pulse is very small at this instant, i.e., Ta > Th, 
and more specifically r > Tn. A condition for diabaticity 
is derivable imposing that the adiabaticity parameter at 
^2 is aad(i2) ^ 1, thereby finding a bound on the time- 
lag T — t2 — tl which also depends on the form of the 
pulse envelopes. 
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